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Abstract 

We study the electronic structure of small lithium clusters Li„ (n = 2 ~ 8) using the electronic 
stress tensor. We find that the three eigenvalues of the electronic stress tensor of the Li clusters 
are negative and degenerate, just like the stress tensor of liquid. This leads us to propose that we 
may characterize a metallic bond in terms of the electronic stress tensor. Our proposal is that in 
addition to the negativity of the three eigenvalues of the electronic stress tensor, their degeneracy 
characterizes some aspects of the metallic nature of chemical bonding. To quantify the degree of 
degeneracy, we use the differential eigenvalues of the electronic stress tensor. By comparing the 
Li clusters and hydrocarbon molecules, we show that the sign of the largest eigenvalue and the 
differential eigenvalues could be useful indices to evaluate the metallicity or covalency of a chemical 
bond. 
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I. INTRODUCTION 



The stress tensors in general are widely used for description of internal forces of matter 
in various fields of science such as mechanical engineering and material science. In quantum 
systems as well, the stress tensors have been investigated for many years, including one of the 
earliest quantum mechanics papers [lH23]. As for the stress tensors in quantum mechanics 
context, we can find several different definitions and applications in the literature. For 
example, Ref. [6] and followers focus on the stress tensor which is associated with forces on 
nuclei. In contrast, the one we consider in this paper is the electronic stress tensor, which is 
associated with effects caused by internal forces acting on electrons in molecules, following 
Ref. [I2j. This electronic stress tensor has been used to investigate chemical bonds and 
reactions and many interesting properties have been discovered [121 [13 [El EDI [2H432] . 

We briefly review our past works on the electronic stress tensor which are closely related 
to this paper. In Ref. [15] , it has been proposed that a covalent bond can be described by the 
eigenvalues and eigenvectors of the electronic stress tensor. In detail, the bonding region with 
covalency can be characterized and visualized by the "spindle structure" , where the largest 
eigenvalue of the electronic stress tensor is positive and the corresponding eigenvectors form 
a bundle of flow lines that connects nuclei. In Ref. [24] , the eigenvalues of the electronic stress 
tensor of various homonuclear diatomic molecules have been investigated and all metals and 
metalloids molecules studied there exhibit the negative largest eigenvalue at the midpoint 
of the nuclei. Ref. [21] has mentioned that this negativity of the largest eigenvalue may 
be connected to the metallic nature of bonding but, at the same time, it has been known 
that the negative largest eigenvalue appears between a pair of atomic nuclei with very short 
distance such as the C-C bond in C2H2 US [2l|. In Refs. [30] and [32], the AI4 and Pd 
clusters have been investigated and we have shown "pseudo-spindle structure" [32] between 
metal nuclei. The pseudo-spindle structure is a negative eigenvalue version of the spindle 
structure and, more precisely, it is the region between two atoms where the largest eigenvalue 
of the electronic stress tensor is negative and corresponding eigenvectors forming a pattern 
which connects them. However, again, the pseudo-spindle structure is also seen in C2H2 
[T6l l3T] . Since the pseudo-spindle structure is seen in many bonds between metal nuclei, 
it is tempting to associate it with a metallic bond, but it seems that the pseudo-spindle 
structure is not sufficient to capture the entire feature of a metallic bond. 
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In this paper, we would like to argue how we may describe a metallic bond using the 
electronic stress tensor. The conventional wisdom would say that it is not sensible to define 
"a metallic bond" because metal is something defined only for bulk with many electrons and 
ions. Our challenge is capturing some aspects of a metallic bond or metallicity of a chemical 
bond using the electronic stress tensor. As a starting point, we tackle this problem using 
small clusters of lithium, the simplest metal. 

There are many studies on the nature of bonding in the Li clusters in the literature 
using such methods as by generalized valence-bond [33] , electron localization function (ELF) 
[SUES] and quantum theory of atoms in molecules (QTAIM) We would also like to 

see whether our stress tensor approach gives complementary views to these methods. 

The structure of this paper is as follows. Sec. ITTlbriefiy summarizes our analysis method 



of electronic structures using the electronic stress tensor. Sec. |III| shows our results of 
stress tensor analysis on the Li clusters. In Sec. |III A[ we discuss Li2. We also discuss H2 



and LiH for comparison. In Sec. IIIB, the Li clusters are analyzed and a method to use 



differential eigenvalues of the stress tensor is proposed. In Sec. Ill C that method is applied 



to hydrocarbon molecules and its effectiveness is tested. Finally, Sec. [IV] is devoted to our 
conclusion. 



II. THEORY AND CALCULATION METHODS 

In this section, we summarize quantities and method we use for our electronic stress 
tensor analysis. Our stress tensor is based on the rigged quantum electrodynamics (QED) 
theory [12]. The rigged QED is the ordinary QED (the theory of photons and electrons) 
equipped with degrees of freedom of atomic nuclei in order to be used for atomic or molecular 
systems. Since this is a theory of quantum fields, there exist quantum operators defined at 
each point in spacetime. By using the equation of motion for operators of various physical 
quantities, many interesting operator relations have been found [121 113 HH EO] . 

The one related to the electronic stress tensor is the equation of motion of the kinetic 
momentum density operator. It turns out that the time derivative of the kinetic momentum 
density operator can be expressed by the sum of the two force operators: the tension density 
operator f^{f) and the Lorentz force density operator L{f). T^{f) can be expressed by 
the divergence of some 3 by 3 tensor operator, which we call the stress tensor density 
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operator r (r*). We obtain the quantum field theory version of the equihbrium equation 
of stress tensor by taking the expectation value of this equation of motion with respect to 
stationary electronic state of an atom or molecule (note that the expectation value of the 
time derivative of any operator with respect to stationary state is zero) as = f^{r) + {L{f)) 
or = J^i^^i'^^'^K^ + (-^'^(^)) where {k,l} = {1,2,3}, (...) denotes the expectation value, 
r'^(f) = (r'^(r)) and r^^^r) = (f'^^'(r)) 

The expectation values for the stress tensor density and tension density can be expressed 

as 
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where {k, /} = {!, 2, 3}, m is the electron mass, ipiir) is the ith natural orbital and z/j is its 
occupation number. Since the stress tensor is an Hermitian matrix, its eigenvalues are real 
and eigenvectors are orthogonal to each other. We denote the eigenvalues by Aj {i = 1, 2, 3) 
with the ordering of A3 > A2 > Ai. 

The concept of "Lagrange point" is proposed in Ref. [21] and we use it as a point which 
can characterize a bond between two atoms. The Lagrange point fi is the point where the 
tension density r'^(r) vanishes: T^^{fL) = 0. We report the eigenvalues of electronic stress 
tensor at this point in the following sections. 

We also analyze electronic structure of molecules using kinetic energy density nT{r), 
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which is defined in Ref. 



Note that our definition of the kinetic energy density is not 
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positive-definite. Using this kinetic energy density, we can divide the whole space into 
three types of region: the electronic drop region with nj'(r) > 0, where classically 
allowed motion of electron is guaranteed and the electron density is amply accumulated; the 
electronic atmosphere region Ra with nT{r) < 0, where the motion of electron is classically 
forbidden and the electron density is dried up; and the electronic interface S with nT{r) = 0, 
the boundary between Rd and Ra, which corresponds to a turning point. The S can give 
a clear image of the intrinsic shape of atoms and molecules and is therefore an important 
region in particular. 

The electronic structures used in this paper are obtained by the Gaussian 09 [32] for clus- 
ter models and ABINIT [lOl Hi] for periodic models. We use Molecular Regional Density 
Functional Theory (MRDFT) package |12] to compute the quantities, Eqs. ([T|-([3]), intro- 
duced in this section. Some part of the visualization is done using PyMOL Molecular Viewer 
program |l3j. 

III. RESULTS AND DISCUSSION 
A. Li2 

In this section, we analyze the smallest Li cluster, Li2, by the electronic stress tensor. 
For comparison, we also show the results of stress tensor analysis of H2 and LiH. The 
computation is performed by the coupled-cluster single-double (CCSD) method using the 
6-311++G** basis set 

Before we discuss Li2, we review how the chemical bond of the hydrogen molecule is 
expressed by the electronic stress tensor (Eq. ([T])). In Fig.[l](b), we plot the largest eigenvalue 
of the stress tensor and corresponding eigenvector on the plane including the internuclear 
axis. If the largest eigenvalue of the stress tensor is positive, it is called the tensile stress, 
and if negative, it is called the compressive stress. Namely, the sign of the largest eigenvalue 
tells whether electrons at a certain point in space exhibit tensile stress or compressive stress. 
The difference between tensile stress and compressive stress are described as follows. If 
we consider a fictitious plane at some point in atomic or molecular system, there are two 
ways that a part on one side of the plane acts on the other part of the plane. One is when 
electrons on one side of the plane is "pulled up" by the electron field on the other side. The 
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other is when electrons on one side is "pushed back" by the electron field on the other side. 
The former corresponds to the tensile stress and the latter to the compressive stress. The 
direction of the eigenvector corresponds to the direction normal to the plane. In Fig. [l] (b), 
we can see that the region with positive eigenvalue spreads between the H atoms and, in 
that region, the eigenvectors form a bundle of fiow lines that connects the H nuclei. Such 
a region is called "spindle structure" [15] and it is clearly seen in the figure. This means 
that electrons close to one H nucleus are pulled up by the electron field close to another H 
nucleus. Such a structure can be associated with the formation of the tight bonding and has 
a suitable directionality for a bond between two H nuclei. Thus, it is proposed that such 
positive eigenvalue region is the manifestation of the strong covalent bond [15]. 

We now turn to Li2, whose largest eigenvalue of the stress tensor and corresponding 
eigenvector is shown in Fig. [2] (b) in the same manner as H2. The striking difference is 
that there is a region with the negative largest eigenvalue between Li nuclei and there is 
no spindle structure. We do not even see a "pseudo-spindle structure" [32], the negative 
eigenvalue version of the spindle structure {i.e. the region with the eigenvectors forming 
a bundle of flow lines that connects two nuclei but with negative eigenvalues), which is 
frequently seen between metallic atoms [301132] . This means that the electron fleld pushes 
back electrons in the neighboring regions with no solid directionality which can be considered 
as a bond axis. This is similar to the stress tensor of liquid. A liquid has the ability to 
flow. Liquid particles are loosely bound and they can move around one another freely. Such 
properties can be well associated with electrons in metal. In fact, coinci dentally, metallic 
bonding is often described as an "electron sea" and the positively charged metal ions. It 
may be unusual to say Li2 has purely a metallic bond with no trace of covalency but if we 
look at the electronic stress tensor, it is totally different from a covalent bond and certainly 
has a feature which can be called metallic. We will return to this point when we discuss 
larger Li clusters. 

After seeing the difference between H2 and Li2 in terms of the electronic stress tensor, it 
may be interesting to see what the stress tensor of LiH look like. This is shown in Fig. |3] (b). 
It has both positive and negative regions between H and Li nuclei. The region close to H/Li 
has positive/negative eigenvalue. In the region very close to the internuclear axis, there 
is some flow of eigenvectors in the direction of connecting the nuclei, but it is somewhat 
different from spindle or pseudo-spindle structures. Therefore, we can say that LiH has a 
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bond with solid directionality like covalent bond of H2 but there is a strong part nearby H 
and a weaker part nearby Li. It is tempting to associate this with an ionic bond, but we 
reserve conclusion until more examples are examined. 

Since we have discussed the analogy between the stress tensors of Li2 and liquid, it is 
important to investigate the other two eigenvalues. We plot all the three eigenvalues of the 
stress tensor for H2, Li2 and LiH along the internuclear axes in Fig. |4j The bond length 
and the eigenvalues at the Lagrange point are summarized in Table |T] Let us first look at 
H2. As shown in Fig. |4] (a), around the midpoint of the two H nuclei, the largest eigenvalue 
(A3) is positive and the smaller eigenvalues (Ai and A2) are negative and degenerate. The 
absolute values of the three eigenvalues are (9(0.1) and so are the differences between A3 
and Ai or A2. Fig. |4] (b) shows the case of Li2 and we see that the three eigenvalues are 
all negative around the midpoint of the two Li nuclei. Two larger eigenvalues (A3 and A2) 
are degenerate and the smallest one (Ai) is smaller than them by only (9(10~^). Therefore, 
compared to the case of H2, the three eigenvalues of Li2 can be regarded as degenerate. Such 
degeneracy in the eigenvalues of stress tensor, suggesting the lack of directionality, solidifies 
the analogy with liquid. Note that, in terms of the stress tensor, liquid is defined by the 
three eigenvalues which are negative and degenerate. This is in stark contrast with the case 
of H2, whose eigenvalue pattern indicates the strong directionality of the bond. The case of 
LiH is shown in Fig. |4] (c). The three eigenvalues are almost degenerate in the region close to 
Li and the region close to H has positive largest eigenvalue and two smaller eigenvalues are 
negative and degenerate. It is different from both Li2 and H2 and it looks like a mixture of 
Li2 and H2. Such a feature could be a hint to consider how an ionic bond may be expressed 
by the electronic stress tensor, but, again, we need more examples. 

So far, we have analyzed the molecules at the equilibrium distances and discussed their 
nature of bonding. We now try to extract more information by changing the internuclear 
distances. In particular, if we compute the stress tensors at very long distances, we can 
study the nature of bonding from the viewpoint of the bond formation. 

The eigenvalues and eigenvectors of the stress tensor for II2, Li2 and LiH at various 
internuclear distances are plotted in Figs. [T| |2] and |3| In the following discussion, the zero 
surface of the kinetic energy density (Sec. [ll| Eq. ([s])), the electronic interface S, plays an 
important role. As is explained in Sec. [Tl| is the boundary between the electronic drop 
region (where the electron density is accumulated) and the electronic atmosphere region 
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Ra (where the electron density is dried up). In the figures, S is drawn by green dashed fines. 
To identify wfiicfi side of S fias positive or negative kinetic energy density {Rd or Ra), it is 
useful to remember Rds exist, naturally, around nuclei. Wfien we look at tfie figures with 
very long internuclear distances, we see S exists around each nucleus, indicating that they 
do not form molecules and suggesting the existence of independent atoms. Moving to the 
figures with shorter internuclear distances, we see that S"s around nucleus get closer, merge 
and eventually form S for the molecules at the equilibrium distances. Incidentally, the state 
when S"s just begin to merge, namely touch each other, is called the "intrinsic electronic 
transition state" [12]. We denote the internuclear distance at this state by r'^^. 

The notable fact is that, when the distance is so long that 5"s are well separated (when 
practically considered as an infinite distance), we see the spindle structure at Ra between 
the two nuclei for all of H2, Li2 and LiH (see Figs, [l] (d), [2] (h) andjs] (d) respectively). This 
means that electrons in i?^) nearby nuclei are pulled up toward Ra through S, resulting 
in the formation of a Lewis electron pair. In other words, all the molecules are covalently 
bonded at the very long distances. In the case of H2, as shown in Fig. [1} the spindle structure 
continues to exist at shorter distances including = 2.20 A and the equilibrium distance. 

In the case of Li2, more interesting things happen. First, we observe that the spindle 
structure exist at the intrinsic electronic transition state {rl^^^ = 5.43 A, Fig.g(g)). Then, 
at 3.31 A, the positive eigenvalue region vanishes and the spindle structure turns into the 
pseudo-spindle structure. (Note that, at slightly longer distance of 3.36 A, the positive 
eigenvalue region vanishes on the internuclear axis while positive areas remain away from 
the axis. Namely, between 3.36 A and 3.31 A, the positive region is not simply connected.) 
That pseudo-spindle structure is broken at 2.78 A, at which three eigenvalues are degenerate 
(Fig. [5] (b)). At shorter distances, we see a negative region with no particular directionality 
such as to connect two nuclei, especially at the equilibrium distance. This is an example of 
what has been argued in Ref. |22]- When there are two distant Li atoms, the tensile stress 
(as is indicated by the spindle structure) pulls up electron in i?/) to the adjacent Ra through 
S which separates them. The consequence is the formation of the long-range Lewis pair of 
electron. As the distant pair of Li atoms comes closer, the spindle structure disappears, 
indicating that the Lewis pair is unbound. This is considered to be a manifestation of 
metallicity. 

In the case of LiH, at long distances, Ra between Li and H is characterized by the spindle 
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structure as mentioned above and, in particular at = 3.38 A. In Figjs] (c), we see the 
contact point of the two interfaces is covered by the spindle structure. When the distance 
gets shorter, although the spindle structure is partly lost, the positive eigenvalue region 
remains around H. The negative eigenvalue region develops in the region close to Li but 
do not cover the whole space unlike the case of Li2. As for LiH, from the viewpoint of the 
electronic stress tensor, we may conclude that it has a bond which is neither covalent nor 
metallic. Whether we can categorize it as ionic or other bonding is an interesting question 
to ask and will be investigated in future. 



B. Lithium clusters 

The optimized structures for Li clusters, Li„ (n = 3 ~ 8) are shown in Fig. [6] In 
the figure, atoms are connected when the Lagrange point (Sec. [ll]) is found between them. 
These structures are obtained by re-optimizing the structures reported in Ref. [S]. We 
performed optimization with the three lowest multiplicities for each cluster and adopted the 
one with the lowest energy, which turned out to be singlet for clusters with even numbers 
of atoms and doublets for those with odd numbers. The computation is performed by the 
CCSD method using the 6-311++G** basis set [l5l HGj. Although computational setups 
are slightly different, the obtained structures are basically consistent with Ref. |44j- As for 
Lie, Ref. [B] has reported the structure with D2h symmetry while we have obtained the one 
with higher symmetry of D^h, which is consistent with Ref. |47j . 

We first examine the largest eigenvalue of the stress tensor and corresponding eigenvector 
on the plane including three atoms in the Li clusters. They are plotted in Figs. [7| and |8j We 
omit equivalent bonds due to the structural symmetry. In these figures, every pair of atoms 
has a Lagrange point in between. We focus on the electronic drop regions {i.e. regions with 
positive kinetic energy density) between two atoms in each panel. The common feature is 
that they all have the largest eigenvalue which is negative. As for the eigenvectors, many 
regions do not show clear pattern such as to connect two atoms, just as the case of Li2. The 
regions with rather clear pseudo-spindle structures are found but not so common (the 1-2 
bond in Lig (Fig. |7| (f)), the 1-2 bond in Liy (Fig. g (a)) and the 1-2 bond in Lig (Fig. |8| (c) 
or (d)). 



From the discussion in Sec. Ill A in order to characterize a chemical bond by the electronic 
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stress tensor, it is important to know the pattern and degree of degeneracy of the three 
eigenvalues in addition to the sign of the largest eigenvalue. The three eigenvalues of stress 
tensor at the Lagrange points are summarized in Table [TTj From the table, we can reconfirm 
that the eigenvalues are all negative. In addition, we see that the three eigenvalues are close 
to each other as is the case of Li2 at least at the Lagrange points. 

In order to quantify the degeneracy pattern more clearly, we propose to use following 
differential eigenvalues. 

Ad32 = A3 — A2, (4) 
^D2i = A2 — Ai. (5) 



Since our convention is A3 > A2 > Ai, A/332 and Xd21 are positive. For example, at the 
Lagrange point, (Ad32, )^D2i) = (0.394, 0) for H2, (0, 2.37 x 10"^) for Lia and (8.52 x lO'^, 0) 
for LiH. In Fig. |9] (a), we plot (Ad32, ^021) at the Lagrange points for the Li clusters. It is 
shown that the Li clusters have very small values of Ad32 and Ad2i, the same order as that 
of A£)2i of Li2 and much smaller than Xd32 of H2. Thus, we can consider three eigenvalues 
of the Li clusters are almost degenerate, meaning the high degree of isotropy (low degree 
of directionality) of the electronic stress tensor. Together with their negativity, the stress 
tensor of the Li clusters are just like that of liquid. As we have mentioned in Sec. Ill A 



electrons in metal are readily associated with particles in liquid. By turning this argument 
around, we may try to consider that negativity and degeneracy of the eigenvalues of the 
electronic stress tensor characterize and quantify some aspects of the metallic nature of 
chemical bonding. 

To support the validity of this idea, we here show some data of very small Na clusters, 
Na„ (n = 2 ~ 4). The electronic structures are obtained by the same computational setups 
as the Li clusters. We performed geometrical optimization starting from the structures 
reported in Refs. |17l SH] . The most stable structures turned out to be singlet for clusters 
with even numbers of atoms and doublets for those with odd numbers. As for the symmetry 
of those structures, Nas and Na4 are found to have and symmetry. The results of 



the electronic stress tensor analysis are summarized in Table III From the data in the table, 
we readily see that the eigenvalues are all negative and almost degenerate just as the case 
of Li clusters. 
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As another way to test this idea, we analyze the electronic stress tensor of periodic 
models for bulk Li and Na. The body-centered cubic structure is adopted for both Li 
and Na with the lattice constants taken to be 3.491 A and 4.225 A respectively. We use the 
norm-conserving pseudopotentials of TrouUier- Mart ins type [19] and the generalized-gradient 
approximation method by Perdew-Burke-Ernzerhof [50] for density functional exchange- 
correlation interactions. Kinetic energy cutoff of plane-wave expansion (k-point) is taken 
as 40.0 hartree (2x2x2 k-point set). We compute the electronic stress tensor at the 
midpoint of two nearest neighborhood atoms. We obtain A3 = —0.468 x 10^^ and A2 = 
Ai = -0.716 X 10-3 for Li, and A3 = A2 = -0.253 x 10"^ and Ai = -0.293 x 10^^ for Na. 
As for the differential eigenvalues, (Ad32, Az)2i) = (2.48 x 10^^,0) for Li and (0,3.95 x 10"^) 
for Na. These results show that the negativity and degeneracy of the eigenvalues of the 
electronic stress tensor are found in the bulk systems too, and the degrees of negativity and 
degeneracy are not different from the clusters. 

Here, some comments are in order as for comparison with the studies in the literature. 
As shown just above, from the viewpoint of the electronic stress tensor, the chemical bonds 
in the small Li clusters look similar to the bulk Li. They posses compressive stress (negative 
eigenvalues) with no particular directionality (almost degenerate eigenvalues) like liquid, 
which can be considered as one of the important features of metal. Interesting thing is that 
this is even so with Li2. Past studies like Refs. |33l[3l] show multi-center bonds are ubiquitous 
features of the Li clusters. In particular, Ref. [31] argues that the metallic phase of Li is 
characterized by multi-center bonds and very gentle variation of ELF across space. From 
this viewpoint, since Li2 is not multi-centered by definition, Li2 does not have a common 
feature with bulk Li or the other Li clusters. This point is in stark contrast with the stress 
tensor viewpoint. We note that, in Ref. [36l 137] , where Li2 is studied in the framework of 
QTAIM, Li2 is viewed as a system that two Li atoms are bonded to a pseudo atom (the 
central nonnuclear attractor) and electrons at the central regions are associated with those 
"partially behave as mobile metallic electron" . The viewpoint of the stress tensor is closer to 
one of QTAIM in a sense that Li2 already exhibits some aspects of chemical bonding which 
characterizes bulk Li. 
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C. Comparison with other molecules 



In this section, we analyze hydrocarbon molecules and compare with the Li clusters, 
especially regarding the negativity and degeneracy of the eigenvalues of the electronic stress 
tensor. Ref. pBj has studied hydrocarbon molecules using the electronic stress tensor. We 
use similar set of molecules but not exactly the same one. The complete list of structures of 
our set which consists of 48 molecules is shown in the supplementary material (Fig. SI and 
Table SI [51]). The structures of these hydrocarbon molecules are obtained by optimization 
using the density functional theory (DFT) method with B3LYP [3^ as a functional and 
the 6-311++G** basis set. 

As for the electronic stress tensor analysis, we repeat the same analyses which has been 
carried out for the Li clusters in the previous section. The detailed data table for the 
hydrocarbon molecules like Table |ll] can be found in the supplementary material (Table S2 
[51j). We begin by investigating the largest eigenvalue of stress tensor. We find that most of 
the bonds have positive largest eigenvalues. Negative largest eigenvalues are found in some 
of the triple C-C bonds, as is mentioned in Sec. |T} In our data set, the triple C-C bonds in 
C2HX with X being H, CII3, C2H5, NH2, OH, F, CI and Br have the largest eigenvalues 
which are negative and whose absolute values are 0(10^^) or less (the three eigenvalues for 
these bonds against the bond lengths are plotted in Fig. S2 in the supplementary material 
[SI]). It is known that this negativity of the largest eigenvalue in these triple bonds can 
be attributed to the very short bond and not to metallicity [TU 123 ISI]- The pattern of 
degeneracy of three eigenvalues of the electronic stress tensor is very different from that of 
the Li clusters. In the hydrocarbon molecules, the smaller two eigenvalues are degenerate, 
Ai ~ A2, and A3 has much larger values which are mostly positive. Thus, if we only focus on 
the largest eigenvalue, these triple C-C bonds and the bonds in the Li clusters look similar. 
However, by looking at the other two eigenvalues, we can clearly distinguish them. 

As is done in the previous section, we compute the differential eigenvalues Xd32 and \d2i 
at the Lagrange points in the hydrocarbons. This is plotted in Fig. |9] (b) and we see that 
Xd32 is much larger than Xd2i- It means that the smaller two eigenvalues are relatively 
degenerate and the largest one is much larger than those two. This distribution pattern 
is quite different from that of the Li clusters shown in Fig. |9] (a). Let us summarize our 
results using the differential eigenvalues. The Li clusters are characterized by Ad32 ^ 1 and 
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Ad2i <^ 1 whereas the hydrocarbon molecules Xd32 ^ ^£121 and Xd2i <^ 1- In particular, 
AD32(Li„) ?s A£)2i(Li„) <S A£)32(h/c). This means that the electronic stress tensor between 
atom pairs in the hydrocarbon molecules shows directionality while it does not show much 
directionality in the case of the Li clusters. 

In summary of this section, first of all, we find that most of the hydrocarbon molecules 
are different from the Li clusters in a sense that the largest eigenvalue of the stress tensor 
is positive. Although some of the triple C-C bonds have the largest eigenvalue which is 
negative, they do not look hke those of the Li clusters due to the non-degeneracy of the 
eigenvalues. There is no hydrocarbon molecule which has the electronic stress tensor like 
that of liquid, at least in our hydrocarbon molecule set. Thus, this results arc not inconsistent 
with the idea proposed in the last section that negativity and degeneracy of the eigenvalues 
of the electronic stress tensor may characterize and quantify some aspects of the metallicity 
of chemical bonding. 

IV. CONCLUSION 

In this paper, we have studied the electronic structure of small lithium clusters Li^^ (n — 
2 ~ 8) using the electronic stress tensor. We have found that the degeneracy pattern of 
the three eigenvalues A, (i = 1,2,3) (A3 > A2 > Ai) of the electronic stress tensor at the 
Lagrange points of the Li clusters is very much different from that of the hydrocarbon 
molecules, which are covalently bonded. Namely, the three eigenvalues of the Li clusters 
have almost same values while the hydrocarbon molecules have the largest eigenvalue much 
larger than the second largest eigenvalue, which has similar value to the smallest eigenvalue. 
The former degeneracy pattern indicates that the bonds are not directional while the latter 
indicates the clear directionality of the bonds. The negativity of the largest eigenvalue is 
not unique feature of the Li clusters (some of the C-C triple bonds exhibit the negative 
largest eigenvalue) but this degeneracy is characteristic to Li. We can consider that such 
difference in the degeneracy pattern reflects the metallic nature of bonding of the Li clusters 
and the covalent nature of bonding of the hydrocarbon molecules. Note that the negativity 
and degeneracy of the stress tensor, which is the property of the liquid, is readily associated 
with the tradition that the metal can be viewed as "electron sea" . 

To describe the degeneracy pattern of the eigenvalues in a compact manner, we have 
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proposed to use the differential eigenvalues: Xd32 — A3 — A2 and Xd2i = A2 — Ai. By 
using them, the degeneracy patterns of the eigenvalues can be summarized as AD32(Li„) fa 
Ai:)2i(Li„) <^ AD32(h/c) ~ C(O.l). It is of great interest whether this relation can be extended 
as Ad32 (metallic) ^ Ad21 (metallic) ^ A£)32(covalent) in general. In this paper, wc have only 
examined the Li clusters with limited numbers of atoms. It is important to check this 
relation with larger Li clusters and also with other metal clusters. Then, the role of the 
electronic stress tensor in describing the metallicity would be elucidated and more detailed 
classification among the "metallic bond" would be possible. 
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TABLE I: The bond length and the eigenvalues Aj (i = 1, 2, 3) of the electronic stress tensor at the 
Lagrange point for H2, Li2 and LiH. 



Bond 


length [A] 


A3 


A2 


Ai 


H2 


0.743 


0.175 


-0.219 


-0.219 


Li2 


2.692 


-1.025x10-3 


-1.025x10-3 


-1.262x10-3 


LiH 


L601 


-4.954x10-3 


-1.347x10-2 


-1.347x10-2 



TABLE XL The bond length and the eigenvalues Aj (i = 1,2,3) of the electronic stress tensor at 
the Lagrange point for bonds in the Li clusters. The numbers in the second column correspond to 
those labelled in Fig. |6| 



Cluster 


bond 


length [A] 


Ac,('xl03) 


A9f xl03) 


Ai(xl03) 


Lis 


1-2 


2.782 


-0.900 


-0.993 


-1.054 




1-3 


3.275 


-0.271 


-0.463 


-0.781 


Li4 


1-2 


2.667 


-0.793 


-0.919 


-1.240 




1-3 


3.010 


-0.498 


-0.763 


-0.981 


Lis 


1-2 


2.929 


-0.516 


-0.708 


-1.152 




1-3 


2.927 


-0.519 


-0.715 


-1.163 




1-4 


3.085 


-0.494 


-0.746 


-0.836 




2-3 


2.618 


-1.139 


-1.157 


-1.348 




2-4 


3.028 


-0.494 


-0.789 


-1.016 




3-4 


3.027 


-0.494 


-0.790 


-1.016 


Lie 


1-2 


2.640 


-1.004 


-1.410 


-1.410 




1-3 


2.849 


-0.600 


-0.924 


-1.198 




1-2 


2.820 


-0.684 


-1.064 


-1.064 




1-3 


2.974 


-0.572 


-0.861 


-0.861 




3-4 


3.078 


-0.465 


-0.636 


-0.787 


Lis 


1-2 


2.926 


-0.584 


-0.826 


-0.918 




1-5 


3.040 


-0.474 


-0.719 


-0.814 
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TABLE IIL The bond length and the eigenvalues Aj {i = 1,2,3) of the electronic stress tensor at 
the Lagrange point for bonds in the Na clusters. Since Naa and Na4 respectively have C2v and 
symmetry groups, which are same as the Li cluster case, one may refer to Fig. [6] for the numbers 
in the second column. Note that we do not find the Lagrange point between the atoms 1 and 3 
(with the separation of 4.371 A) of Nas. 



Cluster 


bond 


length [A] 


A3(xl03) 


A2(xl03) 


Ai(xl03) 


Nas 


1-2 


3.162 


-0.532 


-0.532 


-0.752 


Nas 


1-2 


3.316 


-0.430 


-0.444 


-0.475 


Na4 


1-2 


3.232 


-0.418 


-0.453 


-0.514 




1-3 


3.578 


-0.264 


-0.334 


-0.424 
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(a) 0.5 A 



(b) 0.743 A 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

x(A) ^ x(A) 

(c) 2.20 A (d) 3.0 A 




FIG. 1: The largest eigenvalue of the stress tensor (color map) and corresponding eigenvector 
(black rods) of H2 at various internuclear distances: (a) 0.5 A, (b) 0.743 A (equilibrium distance), 
(c) 2.20 A (intrinsic electronic transition state) and (d) 3.0 A. As for the eigenvectors, the projection 
on this plane is plotted. The red solid line denotes a contour where the eigenvalue is zero. The 
green dashed line denotes a contour where the kinetic energy density is zero (electronic interface). 
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(a) 1 .5 A 



(b) 2.69 A 
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FIG. 2: The largest eigenvalue of the stress tensor and corresponding eigenvector of Li2 at various 
internuclear distances (a) 1.5 A, (b) 2.69 A (equilibrium distance), (c) 2.78 A, (d) 3.31 A, (e) 3.36 A, 
(f) 4.0 A, (g) 5.43 A (intrinsic electronic transition state) and (h) 6.0 A. Plotted in the same manner 
as Fig. [ij 
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(a) 1 .0 A 



(b) 1 .60 A 
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FIG. 3: The largest eigenvalue of the stress tensor and corresponding eigenvector of LiH at various 
internuclear distances: (a) 1.0 A, (b) 1.61 A (equilibrium distance), (c) 3.38 A (intrinsic electronic 
transition state) and (d) 4.0 A. Plotted in the same manner as Fig. [l| 



22 



(a) H2 

0.4 



0.3 
0.2 

tn 0.1 
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tu 
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iU -0.1 
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(b) Li2 

0.001 



(c) LiH 
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x[k] 



-1 1 

x[k] 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 

x|A] 

FIG. 4: Three eigenvalues of the stress tensor for H2 (panel (a)), Li2 (panel (b)) and LiH (panel 
(c)) along the internuclear axis. The red solid lines are for the largest eigenvalue (A3), the green 
dashed lines are for the second largest eigenvalue (A2) and the blue dotted lines are for the smallest 
eigenvalue (Ai). The internuclear distances are set to be their equilibrium distances. The origin of 
the abscissa corresponds to the midpoint of two nuclei. 
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(a) 1 .5 A 



(b) 2.78 A 



0) -0.012 
"(3 



-0.013 
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(d) 6.0 A 




r[Al 



r[Al 



FIG. 5: Three eigenvalues of the stress tensor for Li2 along the internuclear axis at various inter- 
nuclear distances: (a) 1.5 A, (b) 2.78 A, (c) 3.31 A and (d) 6.0 A. The red solid lines are for the 
largest eigenvalue (A3), the green dashed lines are for the second largest eigenvalue (A2) and the 
blue dotted lines are for the smallest eigenvalue (Ai). 



24 



(a) Lis iC2v) (b) Li4 {D2h) 




(c) Lis (Cs) (d) Lie (^W 




(e) UyiDsh) (f) Lis (7-d) 




FIG. 6: Optimized structures of Li clusters: (a) Lis {C2v), (b) Li^ {D2h), (c) Lis iCs), (d) Lig 
{Dih)-, (e) Liy {D^h) and (f) Lig (T^). The bonds are drawn at which the Lagrange points are 
found. 
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FIG. 7: The largest eigenvalue of the stress tensor (color map) and corresponding eigenvector for 
Li clusters Li„ (n = 3 ~ 6). Plotted in similar manner to Fig. [T] The filled circles denote the 
positions of the atoms and the numbers in parentheses correspond to those labelled in Fig. |6j 
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FIG. 8: The largest eigenvalue of the stress tensor (color map) and corresponding eigenvector for 
Li clusters Li„ (n = 7 ~ 8). The optimized structures are used. Plotted in similar manner to 
Fig. [3 
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FIG. 9: The relation between the differential eigenvalues, Xd32 = A3 — A2 and Xd21 = A2 — Ai, 
at the Lagrange points in the Li clusters (panel (a)) and hydrocarbon molecules (panel (b)). Note 
that the H2 molecule has Xd32 = 0.394 and Xd21 = 0. 
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